Off equilibrium properties of vortex creep in superconductors 
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We study a model for the dynamics of vortices in type II 
superconductors. In particular, we discuss glassy "off equi- 
librium" properties and "aging" in magnetic creep. At low 
temperatures a crossover point is found, T g , where relaxation 
times seem to diverge a la Vogel-Tamman-Fulcher. Magnetic 
creep changes by crossing T 9 : above T g power law creep is 
found asymptotically followed by stretched exponential satu- 
ration; below T g the creep is logarithmic and vortex motion 
strongly subdiffusive. In this region violation of time trans- 
lation invariance is found along with important dynamical 
scaling properties. 

Despite apparent structural differences, the dynamics 
of vortices in superconductors has many common 
features with glass formers such as supercooled liquids, 
spin glasses or polymers J^Q. In superconductors even 
simple quantities such as the sample magnetisation, M, 
exhibit strong "glassy" behaviours and a variety of "his- 
tory dependent" effects has been discussed Com- 
mon examples are hysteretic magnetisation loops which 
at low temperature exhibit "memory" effects: the value 
of M for a given applied field, H, depends on the sweep 
rate of H, i.e., on the history of the sample |fj. "Aging" 
has been also recently observed in magnetic relaxation 
in type II superconductors jfj). Interestingly, these phe- 
nomena are very similar to those found in other glassy 
systems such as supercooled liquids. For instance, after 
a quench, "cooling rate" dependences and "aging" effects 
are usually recorded in particle density measurements || . 

Those similarities are striking when one considers the 
differences which exist at a microscopic level between 
these systems. One is apparent: the interaction ranges 
of vortices in superconductors are typically much larger 
than those of molecules in glass formers. In this perspec- 
tive our aim is to discuss the above connections and relate 
off equilibrium dynamics of vortex matter to the dynam- 
ics of glass former liquids, where relevant progresses have 
been recently accomplished [||J|]. In all these systems, the 
glassy features are due to some general mechanisms and 
unrelated to specific material parameters. Thus, a sta- 
tistical mechanics approach, disregarding sample details, 
seems appropriate to describe them. For instance, many 
of the features of glassy behaviours in fragile glasses, such 
as colloidal suspensions, are well described in the context 
of spin glass models . 

Extensive work has been done on the equilibrium or 
close to equilibrium properties of the so-called vortex 
glass phase proposed by M.P.A. Fisher (see ref.s in 
Here we discuss, in particular, the dynamics of the sys- 
tem when far from equilibrium, a typical situation due to 
the enormous relaxation times found in the low temper- 



atures region. We consider a simple tractable model for 
vortices which can be well understood in the context 
of standard statistical mechanics of disordered systems 
0. Such a model was already shown to successfully de- 
pict a unified picture of creep and transport phenomena 
in vortex physics, ranging from magnetisation loops with 
"anomalous" second peak, logarithmic relaxation at van- 
ishing T , Bean profiles, to history dependent behaviours 
in vortex flow and I-V characteristics, to the reentrant 
nature of the equilibrium phase diagram |7j . 

As mentioned, here we study in details the system off- 
equilibrium creep dynamics. The model gives the follow- 
ing scenario for vortex dynamics. At low temperatures a 
crossover point is found, T g , where relaxation times be- 
come longer than the typical observation scale pi] ]. Ac- 
tually, they seem to diverge at a lower temperature, T c , 
where an "ideal" thermodynamic glass transition can be 
located (in 2d T c is numerically zero). Magnetic creep 
changes when T g is crossed: above T g power law creep 
is found asymptotically followed by stretched exponen- 
tials; below T g the creep is logarithmic and shows strong 
form of "aging" , such as violation of time translation in- 
variance with apparent dynamical scaling properties. Mi- 
croscopically, T g corresponds to a drastic change in vor- 
tex motion which becomes strongly subdiffusive. We also 
outlines the connections of vortex off equilibrium dynam- 
ics with other glassy systems. 

We study a Restricted Occupancy Model (ROM) 0, 
a statistical mechanics model which in the limit of zero 
temperature and infinite upper critical field reduces to 
a cellular automaton introduced in JlO| to study vor- 
tex avalanches. The model is a coarse grained version 
JlO|,[?| of a system of straight vortex lines: a set of re- 
pulsive particles diffusing in a disordered pinning envi- 
ronment. The coarse graining length, 1$, is taken to be 
of the order of the natural screening length of the prob- 
lem (which in our case is the London penetration length 
A). Thus, the coarse graining technique (l2) reduces the 
original interaction potential to an effective short ranged 
one. The price to pay is a loss of information on scales 
smaller than Iq and the necessity to allow multiple oc- 
cupancy on the lattice sites of the coarse grained sys- 
tem. The ROM model is defined by the Hamiltonian 

@ : ft = \ niAijUj - 1 Yh A ™ n i ~ Y^i K n ^ together 
with the constraint < rij < N C 2 on the integer occu- 
pancy variable rii representing the number of vortices of 
the coarse grained cell corresponding to site i (bounded 
by the upper critical field). 

The first term in TL represents the repulsion between 
particles E|. For the above considerations on vortex in- 
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teractions, in the present coarse grained representation 
(where 1$ ~ A) a finite range potential must be consid- 
ered. For simplicity we choose An — Aq; Aij — A\ if 
i and j are nearest neighbours; Aij — for all others 
couples of sites. The second term in 7i concerns the par- 
ticle self-interaction energy and the third one is a pinning 
potential, with a given distribution P(A P ) 1 acting on a 
fraction p of lattice sites (below p = 1/2). For simplicity 
we choose: P(Ap) = (1 - p)5(A*>) + p 5(Ap - A p ). A Q 
sets the energy scale and we choose A — 1.0; Aq = 0.3; 
N c2 = 27; k* ee Ax/Aq € [0,0.3] (below we typically dis- 
cuss data for k* = 0.28). The relation of the parameters 
of the model to material parameters is shown elsewhere 
0, here we only recall that k* is an increasing function 
of the Ginzburg-Landau k and N C 2 oc H C 2- 

Through its surfaces, the system is in contact with 
an external reservoir of "particles" , which schematically 
corresponds to the applied field present in magnetic ex- 
periments on superconductors. Particles are introduced 
and escape the system only through the reservoir, which, 
by definition, has a given density N ext . We performed 
Monte Carlo Kawasaki (MC) simulations (T^| on a square 
lattices of linear size L (L £ {8, 128}) in presence of 
a thermal bath of temperature T. Our numerical statis- 
tical averages run, according to system size, from 128 to 
512 pinning realizations. 

The magnetisation in the system at time t Jl4| is de- 
fined as: M(t) = N in (t)-N ext (t), where N in = £J n % /L 2 
is the overall density inside the system and N ext is the 
applied field. As in usual zero field cooled experiments, 
we record the isothermal relaxation of M(t) after ramp- 
ing at a given "sweep rate" , 7, the field from zero up to 
a given working value N ext . 

The presence, at low temperature, of sweep rate de- 
pendent hysteretic cycles, slowly relaxing magnetisation, 
and similar effects 0, indicate that our system, on the 
observed time scales, can be far from equilibrium. Off- 
equilibrium behaviour is appropriately described by two 
times correlation functions. Thus, along with M(t), we 
also record the magnetic correlation function (t > t w ): 

C(t,t w ) = ([M(t) - M(t w )} 2 ) . (1) 

Since the general dynamical features of C and M are 
very similar, for simplicity we mainly discuss C which 
gives richer information. 

At not too low temperatures, for instance at T = 1.0, 
the system creep is characterised by finite relaxation 
times, but the dynamics is already highly non trivial. 
The two times correlator C(t, t w ) is plotted in Fig.|l| after 
a ramp with 7 = 10~ 3 at T = 1.0. At long times, C(t, t w ) 
is well fitted by the so called Kohlrausch- Williams- Watts 
(KWW) law (see Fig .[!]): 

C(i,^)^C 00 {l-e-[( t -^)M 3 } (2) 

The KWW decay is observed in the asymptotic relax- 
ation of superconductors as well as in glass formers (their 



so called a-relaxation [||) above the glass transition ||. 
The time scale r and the Kohlrausch-exponent j3 depend 
on the temperature T (see Fig.||) and on the overall field 
N ext . The pre-asymptotic dynamics (i.e., t,t w << r) is 
also interesting and characterised by various regimes. In 
particular, for not too short times, a power law is ob- 
served over several decades (see Fig.|l|): 

C(t,t w )^C (^^-J (3) 

The exponent a is almost independent of N ext , a ~ 1.7, 
except at very small or high fields. Notice that the r in 
eq.(|^) is the same as in eq.(0), but the exponents a and 
(3 are numerically different (see Fig.||). 

Interestingly, the above finding in the model of power 
laws followed by KWW relaxations is typically observed, 
for not too low T, in superconductors |2|. This behaviour 
is also typical of supercooled liquids, where the power law 
regime is called the /^-relaxation ||. 

At T = 1.0 or, generally, at not too low temperatures, 
no "aging" is seen: C(t,t w ) is a function of t — t w . This 
is clearly shown in Fig.[l], where we plot C(t, t w )/C 00 for 
several different values of N ext as a function of the scaling 
variable (t — t w )/r: all the data (for all t w and N ext ) fall 
on the same master function (which is more general than 
the above KWW fit). 

The scenario described for T = 1.0 is found in a broad 
region at low temperatures. However, around T = 0.5, 
P drastically decreases and, at the same time, a steep 
increase of r is found (see Fig.||). For instance, at 
Next = 10 for temperatures below T g ~ 0.25, the char- 
acteristic time gets larger than our recording window. 
Thus, below T g (N ext ) the system definitely loses contact 
with equilibrium during our observation and, as shown 
in details below, typical glassy phenomena, such as "ag- 
ing" , are observed. The crossover temperature T g is itself 
a function of 7. It has a physical meaning similar to the 
so called phenomenological definition of the glass transi- 
tion point in supercooled liquids. Exploiting this analogy 
we will call this temperature the glass temperature de- 
spite the fact that in glassy systems it is loosely defined 
S. The presence of an underlying "ideal" glass transi- 
tion point, T c (N ex t), is a non trivial possibility which, 
in many cases (as supercooled liquids) is still under de- 
bate Hl5l. T c is often located by some fit of raw r data 
collected in the high T regime || . 

In our model, in the region where r has a steep increase 
a Vogel-Tamman-Fulcher law (VTF) fits the data 0] 
(see inset of Fig.||): 

T = ToCxp (r^k) (4) 

For example, at N ext — 10, the characteristic time tq is 
very large, r = 1.6 -10 3 , and the characteristic activation 
energy, Eq , is ten times larger than T c : Eq = 1 . 1 and T c = 
0.1, thus the above fit is consistent with an Arrhcnius low 
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(i.e., with T c = 0). The presence of a strong increase of r 
close to a power law or a VTF law is again a mark of the 
apparent similitude with glassy features of supercooled 
liquids and glasses ||. 

Since below T g relaxation times are huge, one might 
expect that the motion of the particles essentially stops, 
apart from their vibration inside cages of other vortices. 
Instead, as shown below, the off equilibrium dynamics 
has remarkably rich properties. In particular, the prop- 
erties of the system slowly evolve with time. These phe- 
nomena are usually summarised by claiming that the sys- 
tems is "aging" . They occur in many different systems 
ranging from polymers, to supercooled liquids j^], spin 
glasses Q or granular media p7[ , and their origin and 
apparent universality are still an open problem @,^,^| . 

In the inset of Fig.|^ we show that C(t,t w ), at T = 0.1, 
exhibits strong "aging" . Notice that C depends on both 
times t and t w : the system evolution is slower the older is 
its "age" t w (here, as in eq. ([[]), t w is the time elapsed from 
the sample preparation at the working field and temper- 
ature). Note the contrast with the case at T = 1.0 of 
Fig.|l|. Such a dynamical "stiffening" is typical of glass 
formers S. It is important to stress that the presence of 
slowly relaxing quantities, such as M(t), must not lead 
to conclude that the system is close to equilibrium jjj. 
In the entire low T region (T < T g ), after a short ini- 
tial power law behaviour, C(t, t w ) can be well fitted by a 
generalisation of a known interpolation formula, often ex- 
perimentally used jl| , which now depends on the waiting 
time, t w : 
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U c \ t w + to 
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We found that to take \i ~ 1 is consistent with our data. 
In the above fit U c / ' \xT only depends on N ext and, in- 
terestingly, to is approximately a linear function of t w : 
to oc t w +tQ, where t^ is a constant. Very interesting is the 
presence of scaling properties of purely dynamical origin 
in the off-equilibrium relaxation. This is shown in Fig|], 
where data for different fields, N ext , and different wait- 
ing times, t w , are rescaled on a single master function. 
The above results imply that for times large enough (but 
smaller than the equilibration time), C(t, t w ) is a univer- 
sal function of the ratio t/t w : C(t,t w ) ~ S(t/t w ). Such a 
behaviour (called "simple aging" p|) is in agreement with 
a general scenario of off-equilibrium dynamics (see Ref. 
[ HI) and has strong analogies with other glassy systems 
j 3|^[ |l7|Jl9[| . Experimental data on C(t,t w ) do not exist 
in vortex matter, but would be extremely important to 
shed light on the real nature of dynamical phenomena in 
superconductors. 

Finally, we stress that in vortex physics transitions 
from low temperature logarithmic to higher temperature 
power law creep are usually experimentally found (see 
references in M ) and give a way to approximately locate 
the position of T g in real samples. 



In the above scenario, microscopic quantities concern- 
ing the internal rearrangement of the system, such as the 
vortex mean square displacement, R 2 {t), can be insight- 
ful. Since vortices can enter and exit the sample, R 2 (t) 
must be properly defined. To this aim, we also made a 
different kind of computer runs where after ramping the 
field from zero to N ex t we close the system (i.e., remove 
the reservoir). Since now the number of particles is fixed 
]l8[ , their positions, fi(t), are well defined at each time 
step and we record: R 2 (t) = j^(J2i(^i(^) ~^(0)) 2 )' were 
N is the total number of present particles. 

R 2 (t) is plotted in Fig.fffor N ext = 10. At high enough 
T, R 2 (t) is linear in t, but at lower temperatures it shows 
a pronounced bending. Finally, below T g , the process at 
long times becomes sub-diffusive: 



R 2 (t)~t v 



(0) 



with v << 1. From this point of view, T g is the location 
of a sort of structural arrest of the system, where particle 
displacement is dramatically suppressed. Interestingly, a 
very similar scenario has been recorded in real supercon- 
ducting samples: for instance in Ref. |Q it was clearly 
shown that vortices are definitely mobile in the low tem- 
perature phase and only "freeze" below a characteristic 
field dependent temperature. 

Summarising, the schematic coarse grained model we 
considered shows magnetic properties very close to those 
experimentally found in superconductors Q. Here, 
in particular, we have focused on the off equilibrium 
creep dynamics and an interesting scenario emerges. At 
low temperatures the system relaxation time r(N extl T) 
grows enormously with an "ideal" divergence a la VTF at 
some T c (N ext ) (interestingly, molecular dynamics simula- 
tions of London-Langevin models seem to confirm these 
results J?|]2l|]). In fact, below a crossover temperature, 
T g > T c , vortex displacement becomes sub-diffusive, sig- 
nalling the presence of a structural arrest. It is impossible 
to equilibrate the system on the observation time scale. 
At low temperatures creep is logarithmic and becomes a 
power law above the crossover T g (N ext ). 

In typical experiments or computer simulations, mag- 
netic properties are measured after ramping the exter- 
nal field at a given rate 7. Whenever 7 is much larger 
than the inverse of the characteristic relaxation time, 
i~(N ext , T), the system is naturally driven off equilibrium, 
simply because it is unable to follow the drive, and "hys- 
teresis", "memory" effects, along with dependences on 
the sweep rate, occur. In the present framework, the 
strict correspondence between the dynamics of vortices 
and other glass formers can be also rationalised. 

Work partially supported by EPSRC, 
INFM-PRA(HOP)/PCI. 
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FIG. 1. Main frame The two time magnetic correlation 
function C(t,t w ) recorded at T — 1.0 for the shown values of 
the external field and waiting times t w . All data are collapsed 
on the same curve when plotted as a function of (t — t w )/r. 
Here r(N ex t,T) is the characteristic creep time. No "aging" 
is present. The bold continuous line is a fit to the KWW 
function of the asymptotic region, and the dotted line is a 
power law fit. Inset The same data as above for N ext = 16 
as a function of t — t w . 
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FIG. 2. The parameters of the Kohlrausch- Williams- Watts 
(KWW) asymptotic relaxation of the magnetisation corre- 
lation, C, as a function of the temperature T, recorded at 
N ex t = 10. Main frame The equilibration time r enor- 
mously grows by decreasing the temperature T. Below the 
crossover temperature T g ~ 0.25, the system relaxation times 
are larger than the observation window. Inset left The 
KWW exponent f3 as a function of T. Inset right Close 
to T g , t plotted as a function of 1/T approximately shows a 
Vogel-Tamman-Fulcher behaviour (see eq.(^)). 



FIG. 3. Inset Logarithmic time relaxation at T — 0.1 
(i.e., below T g ) of the two-times correlation function, C(t,t w ), 
recorded at N ex t = 16 (a value close to the 2nd peak in 
magnetic loops). Main Frame Off equilibrium dynami- 
cal scaling. Superimposed on the same master function, 
1 — (1 + x)^ 1 ^ (m _1 ~ 1), are relaxation data of C(t,t w ) 
recorded for N ext = 4, 10, 16 for each of the shown t w . The 
asymptotic scaling is C(t,t w ) ~ S(t/t w ). 
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FIG. 4. The vortex mean square displacement R 2 (t) at 
N ex t = 10 for several temperatures. Below T a ~ 0.25, R 2 {t) 
is strongly subdiffusive: R 2 (t) ~ t" with v < 1. Straight lines 
are guides for the eye. 
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